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Magnetic excitations in the geometric frustrated multiferroic CuCr02 



O 

(N 

< 



C/3 



M. FrontzckiQ J- T. Haraldsen^-^'-*^ A. Podlcsnyak\ M. Matsuda\ A. D. Christianson\ R. S. 
Fishman^, A. S. Scfat^, Y. Qin^'^, J. R. D. Coplcy^ S. Barilo^ S. V. Shiryaev^ and G. Ehlcrsi 
^Neutron Scattering Science Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
^ Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
'^Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
^ NIST Center for Neutron Research, National Institute of Standards and Technology, Caithersburg, Maryland 20899, USA 

^Department of Materials Science and Engineering, 
University of Maryland, College Park, Maryland, 20742, USA 
^Institute of Solid State and Semiconductor Physics, Minsk 220 072, Belarus 
(Dated: September 1, 2011) 

In this paper detailed neutron scattering measurements of the magnetic excitation spectrum of 
CuCr02 in the ordered state below Tni = 24.2 K are presented. The spectra are analyzed using 
a model Hamiltonian which includes intralayer-exchange up to the next-next-nearest neighbor and 
interlayer-exchange. We obtain a definite parameter set and show that exchange interaction terms 
beyond the next-nearest neighbor are important to describe the inelastic excitation spectrum. The 
magnetic ground state structure generated with our parameter set is in agreement with the structure 
proposed for CuCr02 from the results of single crystal diffraction experiments previously published. 
We argue that the role of the interlayer exchange is crucial to understand the incommensurability 
of the magnetic structure as well as the spin-charge coupling mechanism. 

PACS numbers: 75.25-l-z, 75.30.Ds, 75.47.Lx, 75.85-|-t 
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14 I. INTRODUCTION 

15 Compounds which exhibit both an ordered magnetic 

16 phase and a ferroelectric phase are termed multiferroics. 

17 Especially the multiferroics where the electric polariza- 
la tion can be controlled with a magnetic field and vice versa 

19 arc of continuing interest due to the potential applica- 

20 tions. The most promising candidates for such control- 

21 lable multiferroic have been found among the materials 

22 with inherent geometric magnetic frustratiouii 

23 Different mechanisms leading to spin-charge cou- 

24 phng that have been discussed in the literature include 

25 the magneto-elastic effect;^ the 'inverse' Dzyaloshinskii- 

26 Moriya interaction,—!^ and electric dipolc induction 

27 through hybridization of p — d orbitals as originally pro- 

28 posed by Arimai^ Spin-charge coupling due to magne- 

29 tostriction can occur in coUinear commensurate mag- 

30 netic structures as for instance observed in i?Mn205, 

31 where i? is a rare earth metali^ If magnetic order 

32 with non-zero chirality exists, which may be commcn- 

33 surate or incommensurate with the lattice, the inverse 

34 Dzyaloshinskii-Moriya (DM) interaction induces (by in- 

35 version symmetry breaking) an electric polarization com- 

36 ponent perpendicular to the spiral axis and the propaga- 

37 tion vectori^ Systems in which this situation is realized 

38 include TbMnOa^-i^ MnW04r^"-i^ RbFe(Mo04)2r^'^ 

39 LiCu202,i^--2i and Ni3V208.^"-21 Spin-charge coupling 

40 through Arima's mechanism requires a proper-screw 

41 magnetic structure where the vector of the polarization is 

42 parallel to the screw axis and to the propagation vector, 

43 CuFe02 is the most prominent examplci^^"— 

44 In this article, we report a detailed analysis of the spin 

45 dynamics of the multiferroic system CuCr02 which has 



46 already been studied using a variety of techniques such as 

47 polarization in applied magnetic and electric fieldsf^i^ 

48 electron spin resonance (ESR), — x-ray emission spec- 

49 troscopy, (XES) '^^'^^ single crystal x-ray diffraction;^ 

50 neutron diffraction^^— and inelastic neutron scatter- 

51 ing4^i^ This system is isostructural to CuFe02 and a 

52 detailed comparison of the two systems is instructive. 

53 In contrast to CuFe02 which becomes multiferroic in 

54 an applied magnetic field — or through doping the Fe- 

55 site with Al^ Ga^ or Rb^, CuCr02 enters the multi- 

56 ferroic state in zero field with the magnetic transition. 

57 In both compounds the magnetic structure in the mul- 

58 tiferroic phase is an incommensurate proper-screw mag- 

59 netic structure. However, the propagation vector found 

60 for CuCr02 with t = (t,t,0) and r = 0.3298(1) is very 

61 close to the commensurate value. Unlike the propagation 

62 vector of CuFe02 which in comparison is very different, 

63 T = (r, r, 3/2) with r = 0.207^ 



64 II. EXPERIMENTAL 

65 A detailed account of the sample preparation was 

66 given previously!^ The trigonal crystal structure (space 

67 group RZm) with lattice parameters a = 2.97 A and 

68 c= 17.110 A was confirmed by x-ray powder analysis of 

69 crushed crystals. Further characterization with respect 

70 to their magnetic properties was done using a SQUID- 

71 magnetometer. The obtained susceptibility curves are 

72 similar to data published previouslyi^i^'^i^ Identiiy- 

73 ing the same characteristic points in the susceptibil- 

74 ity data as Kimura et al.— the same two characteris- 

75 tic phase transition temperatures, Tni = 24.2 K and 

76 T^2 = 23.6 K, were obtained for our samples. The Curie- 
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77 Weiss fit between 148 K and 287 K of the inverse suscep- 

78 tibility gave an asymptotic paramagnetic Curie temper- 

79 ature of -200(1) K and an effective moment of 3.88(1) hb 

80 per Cr^+ ion. Measurements of the magnetization mea- 

81 sured along three orthogonal directions, [110], [110] and 

82 [001], are shown in Fig. [T] below. A phase transition at 

83 ffflop ''^ 5.3 T can be seen in these data (the value is 

84 determined from the center of gravity of the peak in the 

85 derivative), in agreement with earlier reports At this 

86 phase transition the electrical polarization is flopped^ in 

87 conjunction with a reorientation of the ordered magnetic 

88 moments 4^ 
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112 performed with an incident energy of 3.53 meV with a 

113 measured resolution of 0.1 meV (FWHM) at the elastic 

114 line. The data obtained on CNCS and DCS have been 

115 reduced using the DAVE software package 



III. THEORY 

The hexagonal symmetry of the CuCr02 lattice pro- 
vides a complex network of possible intra- and inter- 
layer supcrcxchangc pathways^ that are described by 
the Heisenberg Hamiltonian 



H 
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FIG. 1. (Color online) Magnetization measurement along the 
three main crystallographic directions in CuCr02 single crys- 
tals at T = 2 K. The inset shows the derivative of the mag- 
netization with a peak at //flop in the [TlO] direction. 

Ten crystals with a total mass of m ^ 0.6 g were co- 
aligned on an aluminum sheet covering an area of approx. 
20 X 20 mm for inelastic neutron scattering experiments. 
The crystals were platelet like with the c-direction nor- 
mal to the plate surface. The horizontal scattering plane 
was HHL. Experiments were conducted at the Cold 
Neutron Chopper Spectrometer (CNCS) at the Spalla- 
tion Neutron Source in Oak Ridgo^, the HB-1 triple-axis 
spectrometer at the High Flux Isotope Reactor in Oak 
Ridge, and at the Disk Chopper Spectrometer (DCS) at 
the NIST Center for Neutron Research (NCNR),^ 

All experiments used a standard orange cryostat in 
a temperature range from 1.5 to ^ 100 K. The CNCS 
measurements were performed in two settings with dif- 
ferent incident neutron energies, 12.1 meV and 3 mcV, 
respectively. The energy resolution at the elastic line was 
0.4350(6) meV full width at half max. (FWHM) at 12.1 
meV and 0.0649(1) meV FWHM at 3 meV, respectively. 
The HB-1 measurements used constant kf ~ 14.7 meV 
which resulted in an effective energy resolution of 1.84 
meV at 7.5 meV. The collimation was 48-60-60-240 with 
two additional pyrolitic graphite (PG) filters to suppress 
higher order contamination. The DCS measurement was 



121 where S, is the local moment on site i. The superex- 

122 change interactions Jij between sites i and j are antifer- 

123 romagnetic when Jy < 0. An overview of the exchange 

124 paths in respect to the lattice is given in Fig. O The 

125 single-ion anisotropy along the x and z axes is given by 

126 Dx^z, where D > Q produces easy-axis anisotropy and 

127 Z) < produces easy-plane anisotropy, respectively. The 

128 three-dimensional magnetic state is constructed by stack- 

129 ing the two-dimensional configurations ferromagnetically 

130 along the c-axis. 

131 Through an energy minimization of the exchange pa- 

132 rametcrs and anisotropy, the magnetic ground state con- 

133 figuration is determined through a classical approach de- 

134 scribed in Ref. [5^ by defining Sz within any hexagonal 
plane as 



Szm = A 



21+1 cos[t^(2/ + l)a;] 



(2) 



136 where the C21+1 harmonics are produced by the easy axis 

137 anisotropy D^- With Ci set to 1, the amplitude A is 




FIG. 2. (Color online) Considered exchange paths in the 
Heisenberg Hamiltonian. 
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138 obtained from the condition that the maximum value of 

139 [^^(R)! equals S. The perpendicular spin components 

140 Sy are given by 

Sy{R) - V5-5,(R)2 . sgn(sin(T,a;)) . (3) 

141 The ordering wavevector and coefficients C21+1 are 

142 determined by minimizing the energy on a large unit cell 

143 of size ^ 10^ a X a X c, where a is the lattice constant 

144 within a hexagonal plane and c is the separation between 

145 neighboring planes. 

146 Based on this magnetic ground state, the spin dynam- 

147 ics are evaluated using a Holstein-Primakoff transforma- 

148 tion, where the spin operators are given by 5*^^ = 5— a|ai, 

149 Si+ = V^Sai, and Si- — •\/2S'a| (a^ and a| are boson 

150 destruction and creation operators). A rotation of the 

151 local spin operators accounts for the non-coUinearity of 

152 the spins 

153 To determine the spin wave (SW) frequencies ujq, 

154 we solve the equation-of-motion for the vectors vq = 

155 [ciq^, Oq''^, a^Q, flq''^, ...], which may be written in terms 

156 of the 2A^x2iV matrix M(Q) asidvq/dt = -[H2,vq] = 

157 M(Q)vq, where N is the number of spin sites in the unit 

158 cellj^ The SW frequencies are then determined from the 

159 condition Det[M(Q) — wq/] = 0. To assure the local sta- 

160 bility of a magnetic phase, all SW frequencies must be 

161 real and positive and all SW weights must be positive. 

162 The SW intensities or weights are coefficients of the 

163 spin-spin correlation function: 

a/3 

164 where a and /3 are x, y, or A more detailed discussion 

165 of this method is contained in Ref. [s^. Notice that mag- 

166 netic neutron scattering measurements (INS) only detect 

167 components of the spin fluctuations perpendicular to the 

168 wavevector Q. The total intensity /(Q,a;) for an INS 

169 scan at constant Q is given by 

IiQ,Lu)^S{Q,Lu)F^cxp{~{Lu^L,Qf/2S^) , (5) 

170 where 6 is the energy resolution and _Fq is the Cr'^^ mag- 

171 netic form factor. 

172 This approach yields additional information on the 

173 magnetic ground state. The magnetic ground state is 

174 not provided for these systems and must therefore be de- 

175 rived from the energy minimization of the Hamiltonian 

176 possible magnetic structures within the ^ 10'^ a x a x c 

177 cell. Therefore, two energetically degenerate states, for 

178 instance commensurate vs. slightly incommensurate, can 

179 be distinguished. 

180 IV. RESULTS 

181 The inelastic excitation spectrum of CuCr02 in the 

182 HH direction as measured at CNCS with Ei = 12 meV 



183 is shown in the upper panel of Fig. [3] Integration along 

184 the L direction was in the range < L < 5 r. 1. u. 

185 (relative lattice units) which is justified by a rather small 

186 dispersion along this direction. Integration along the per- 

187 pendicular HH direction was within ±0.025 r. 1. u. (cor- 

188 responding to ±2.5 deg. out of the scattering plane). For 

189 comparison the model calculation is shown in the lower 

190 panel. 



Hin (H, H, L) (r. I. u.) 
-1/3 1/3 2/3 
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FIG. 3. (Color online) Upper panel; Magnetic excitation 
spectrum in S(Q,aj) of CuCr02 measured at T = 2 K at 
CNCS. Integration range along L was from to 5 in r. 1. u., 
and along the HH direction ±0.025 r. 1. u.. The intensity 
around _ff = at low energy originates from the halo of the 
primary beam. Lower panel: Spin waves computed from the 
best theoretical model, the modes discussed in the text are 
marked Q, /3. 

191 The low energy mode a originates from the magnetic 

192 Bragg peak in the vicinity of = 1/3 and flattens off 

193 at around 5 meV. It has a cusp like local energy mini- 

194 mum at the magnetic zone boundary at H ~ 1/6. The 

195 intensity of this mode is strongest in the vicinity of the 

196 Bragg peak and falls off towards the zone boundary. This 

197 mode is mainly influenced by the model parameters J2, 

198 J3, Dx and Dz (see above). The minimum of the a mode 

199 a.t H = 1/6 is of considerable interest. It can only be 

200 modeled with the inclusion of an antiferromagnetic next- 

201 next nearest neighbor exchange interaction J3. If J3 is 

202 neglected or ferromagnetic, the excitation would be flat 

203 at i? = 1/6 or would show a local maximum. Analyz- 

204 ing the intensity of the a mode at the zone boundary. 
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205 the measurement shows more intensity at H = 1/2 than 

206 at H = 1/6. In the modehng this leads to a negative 

207 in-plane anisotropy constant (otherwise the intensity 

208 would be higher at H = 1/6). In return, this leads to 

209 a ground state with a proper screw magnetic structure 

210 rather than a cycloid. 

211 The non-zero anisotropy terms and mean that 

212 the a mode must be gapped. The gap is too small to 

213 be unambiguously detected at = 12.1 mcV. However, 

214 with improved energy resolution = 3 meV) a gap 

215 of ^ 0.5 mcV is clearly seen as shown in Fig. 2] Here 

216 the integration along the L-direction is only for a small 

217 range around L = 1. The absolute values of Dx and Dz 

218 are adapted in the theoretical calculations to accurately 

219 model this gap. 

220 An overall weaker and flat /3 mode is observed between 

221 5 and 8 meV. The measurement did not resolve whether 

222 a crossing of the a and /3 mode occurs as suggested by the 

223 calculation, mainly due to insufficient resolution. The (3 

224 mode has a maximum of ~ 7.5 meV at the magnetic zone 

225 boundaries at H — 1/6 and H = 1/2. The energy of the 

226 /? mode at these points is mainly determined by J2 and 

227 to a lesser degree by J3. Kajimoto et al4^ ascribed the /3 

228 mode (referred to as "flat component") to the existence 

229 of an interlayer exchange interaction Jz which is incon- 

230 sistent with our data. In the lower panel of Fig. 3, the 

231 computed spin wave excitation spectrum form the best 

232 theoretical model is shown. The a and 13 mode in this 

233 energy range determine J2 and J3 as well as Ji to which 

234 all parameters are relative. In agreement with data from 

235 the literature r^i^ a survival of magnetic collective dy- 

236 namics up to several times Tn is observed at the position 

237 of the a mode. 



CNCS, r=2K 10 




Hin (H, H, 0.5) (r. I. u.) 



FIG. 4. (Color online) Magnetic excitation spectrum of 
CuCr02 measured at T = 2 K at CNCS with 3 meV incident 
energy. The inset shows a constant-Q cut along the excita- 
tion. Error bars represent ±lcr from counting statistics. 

238 The spin-wave spectrum along the L-direction is 

239 dispersion-less for energies above 0.5 meV as already 

240 mentioned above. However, below the energy gap of 0.5 



241 meV a modulation can be seen Fig. [5] For an energy 

242 transfer of 0.2 meV, the measured intensity along L is 

243 higher at the position of the magnetic Bragg peaks com- 

244 pared to the position between. This intensity pattern can 




1 2 3 4 5 
L in (1/3, 1/3, L) (r. I. u.) 



FIG. 5. Magnetic excitation spectrum in S{Q,ijj) of CuCr02 
measured at T = 2 K at DCS with 3.55 meV incident energy. 
The data is integrated in the HH of 0.32 to 0.34 r. 1. u. from 
the central detector bank. The intensity is color coded in a 
linear scale with the exception of the elastic Bragg peaks with 
two orders of magnitude higher intensity. 

246 

247 be reproduced with the introduction of a ferromagnetic 

248 interlayer coupling Jz ■ The magnitude of the interlayer 

249 exchange is small as is the effect on the excitation spec- 

250 trum. 

251 The data presented so far allow the determination 

252 of the values for the exchange interaction and the 

253 anisotropy terms within the given model. The calcula- 

254 tions replicate satisfactorily the a and /3 excitation modes 

255 as shown in the lower panel of [3] The intensity pattern 

256 of the DCS measurement (Fig. [5]) is modeled with the 

257 small interaction term Jz- The interlayer exchange Jz 

258 also results in the magnetic ground state with the in- 

259 commensurate ordering wavevector — 0.329. With- 

260 out the interlayer exchange the magnetic ground state 

261 would be commensurate. The model Hamiltonian also 

262 reproduces the gap in the excitation spectrum, using the 

263 anisotropy terms, which as a consequence leads to the 

264 splitting of the otherwise degenerated magnetic ground 

265 state. This splitting of the degenerate ground state gives 

266 rise to another excited state at higher energies, with 

267 a spin wave dispersion that mirrors the /S mode from the 

268 ground state but which has an additional gap of 2.2 meV. 

269 The intensity of this mode is weaker than the excitations 

270 from the ground state and cannot be seen in the CNCS 

271 data, likely because, by way of how the (Q,a;) space is 

272 mapped in a time-of-flight measurement with the chosen 

273 settings, only i > 1 is covered at Huj > 8 meV. 

274 Figure [S] shows a contour map of the measurements 

275 taken at HB-1. These are constant- scans with an en- 
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276 ergy difference of 0.5 meV in the range from 1.5 meV 

277 to 15 meV. The measurements are along the {HH2) di- 

278 rection. In this figure, it can be seen that another mode 

279 with nearly the same dispersion exists above the /? mode, 

280 which we identify with the /?' mode resulting from the 

281 calculations. The coarser energy resolution of HB-1 leads 

282 to a partial blur of the /3 and /?' mode. The calculation 

283 yields a gap between both modes of 2.2 meV at the zone 

284 boundary. 



H in (H, H, 2) (r. I. u.) 
1/6 1/3 1/2 




1/6 1/3 1/2 

Hin (H, H, 2) (r. I. u.) 



FIG. 6. (Color online) Upper panel; Contour map from 
constant-i5 scans of CuCr02 measured at T = 2 K at HB- 
1. Lower panel: The corresponding model of the a, j3 and /?' 
excitations. 

285 To summarize the results, the intensity and dispersion 

286 of experimentally observed spin-wave modes in CuCr02 

287 have been modeled with a Hamiltonian that includes at 

288 least six free parameters, which are given in Table HI 



Set 


Ji 


J2 


J3 








Ref. 47 


-2.3 


-0.12 






-0.4* 


0.4* 


This work 


-2.8 


-0.48 


-0.08 


0.02 


-0.59 


0.48 


CuFe02 


-0.23 


-0.12 


-0.16 


-0.06^ 




0.22 



TABLE I. Comparison of the relevant exchange interaction 
and anisotropy parameters from Ref. [13 (*only one value 
was fitted) with this work and the results for CuFe02 from 
Ref. Isrlf^ J^i ). Energies are in meV. 

289 
290 

291 Small discrepancies between calculation and mcasure- 

292 ment suggest the need to include higher order parameters 

293 beyond the ones used here. This is most apparent in the 



294 slight discrepancy of the spin-wave velocities. The ve- 

295 locities depend in a non-trivial way from all interactions 

296 and deviations from the model may indicate the need for 

297 magneto-elastic or bi-quadratic terms. While the addi- 

298 tion of J3 and helps reduce this difference, it is clear 

299 that other interactions may be affecting the system. The 

300 deduction of the parameters in the Hamiltonian has been 

301 based on the approach to incorporate the least necessary 

302 number to describe the excitation spectrum satisfactorily. 

303 In comparison to CuFe02; the nearest neighbor in- 

304 tralayer exchange interaction Ji is one order of magni- 

305 tude stronger in CuCr02, but the interlayer exchange 

306 and the anisotropy parameter Dz are of comparable mag- 

307 nitude.— The different magnetic ground states are ex- 

308 plainable with the different ratio of D/|Ji|. In CuCr02, 

309 where this ratio is small, the proper-screw is the stable 

310 magnetic structure, while in CuFc02 the four-sublattice 

311 coUincar structure is the ground state.— It has been in- 

312 terpreted that the main effect of doping in CuFe02 is the 

313 decrease of anisotropy and through this the proper-screw 

314 magnetic structure can be stabilized as ground state in 

315 the doped compounds)^ Notably is the difference of the 

316 in-plane anisotropy Dx which is absent in CuFe02 where 

317 a Goldstone mode at the incommensurate wavevector is 

318 observed^, but present in CuCr02 as indicated by the 

319 gap of the a mode. Instead of the observed lattice 

320 distortion in the basal plane is relevant to model the ex- 

321 citation spectra in CuFe02i^ 

322 The interlayer exchange in CuFe02 leads to a 10-sub 

323 lattice stacking sequence along the c-direction and can be 

324 modeled with one ferromagnetic and two antiferromag- 

325 netic exchange parametersi^ The interlayer exchange in 

326 CuCr02 seems simpler and can be described with one fer- 

327 romagnetic parameter of similar magnitude. In CuFe02 

328 the interlayer exchange has been the most affected pa- 

329 rameter by doping^ which might explain the difference 

330 between CuCr02 and CuFe02. 

331 The last marked difference to be discussed is the appar- 

332 ent absence of a structural phase transition in CuCr02. 

333 Strain measurements on CuCr02— indicate strong mag- 

334 netoelastic coupling, but apparently insufficient to lead 

335 to a phase transition as in CuFc02. In the latter, it 

336 has been demonstrated that the inclusion of bi-quadratic 

337 terms in the Hamiltonian are relevant in the prediction of 

338 the phase diagrami^ In CuCr02, the bi-quadratic terms 

339 seem less relevant for the understanding of the magnetic 

340 ground state but probably cause the slight discrepancy of 

341 the spin-wave velocities between model and experiment. 



342 V. CONCLUSION 

343 A detailed investigation of the magnetic excitation 

344 spectrum of CuCr02, at low temperatures has been per- 

345 formed using neutron scattering techniques. The exci- 

346 tation spectrum has been used to deduce the relevant 

347 exchange interaction and anisotropy parameters. The 

348 parameter set points to a ground state with an incom- 
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